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Abstract 

In a recent paper, Phys. Rev E 81, 041137 (2010), the author attempts to derive ten neces- 
sary conditions for stability of dissipative fluids and plasmas. Assuming the validity of the local 
equilibrium principle, these criteria have been obtained solely from the first and the second laws 
of thermodynamics. The Onsager reciprocity relations have not been invoked and author's results 
are supposed to be valid independently of the choice of the boundary conditions. In the present 
paper, in agreement with the general theory established by Glansdorff-Prigogine in 1954 and 1970, 
we shall show that there is no variational principle expressing the necessary conditions for stability 
of dissipative systems involving convective effects when the system is out of the Onsager region. In 
particular, we shall prove that the basic equations constituting the starting point of the analysis 
of the author, attempting to derive ten necessary conditions for the stability involving magneto- 
hydrodynamical effects, are incorrect and in contradiction with the laws of the thermodynamics of 
irreversible processes. 
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I. INTRODUCTION 



In the recently published papers fl, Q and Q, the author attempts to derive ten necessary 
stability criteria for dissipative systems involving magneto-convective effects. These criteria 
take the form of - or they are derived from - variational principles and they have been 
obtained by invoking no Onsager symmetry and no detailed model for heat production 
and transport. By integrations of the inequalitiy expressing the stability condition, and 
the balance equations for mass and energy, the author derived a set of constraints on the 
evolution of smooth perturbations relaxing towards the steady-state. These constraints take 
the form of inequalities involving the total time derivative of quantities such as the volume 
of the system and the entropy produced by heating processes. Each inequality takes the 
form d t A < d t B + d t C, where d t stands for the total (i.e., substantial) time derivative and 
A, B and C are volume integrals. 

However, these results are manifestly in contradiction with the laws of the thermodynamics 
of irreversible processes. The mistake in the work published in Refs [l], 01, is a result of 
a misinterpretation of Eqs (1.5) and (1.4), reported in Ref. [lj and Ref. [2], respectively, 
ndeed, these equations are valid only for dissipative systems in absence of convective effects 
Q- We would like to clarify that, contrary to the claims of Ref. [lj], Eqs (1.5) and (1.4) do 
not correspond to the expression of the Universal Criterion of Evolution (UCE) demonstrated 
by Glansdorff and Prigogine in 1954. Eqs (1.5) and (1.4), are derived from the Gibbs 
expression, written in local form, in conjunction with the local equilibrium stability. Hence, 
these equations are a direct consequence only of the second principle of thermodynamics. 
In the case of pure dissipative processes, the time derivative appearing in these equations 
should be understood as the total time derivative (in the substantial sense). However, for 
inhomogeneous systems in presence of convective effects, Eqs (1.5) and (1.4) should be re- 
written by involving only the partial time derivatives, dt, and by no means in terms of total 
time derivatives. In addition, in the presence of convective effects, the inequality for the 
necessary condition of stability should include the extra term — T^ 1 (d t y) 2 Q], where T and 
v are the temperature and the velocity of the matter, respectively. This term has been 
omitted in Eqs (1.5) and (1.4) of Refs lj and 2|. In addition, if q indicates the heat current 
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associated with the heat balance equation of the form 

dT 

pCv-^ = -V ■ q (1) 

with p and c v denoting respectively the mass density and the specific heat at constant 
volume, it becomes obvious that Eq. (2.6) in Ref. cannot be correct. Eq. ([I]) leads to 
the inequality 

where f2 is the volume occupied by the system. However, expression (2.6), deduced in Ref. |2[, 
exhibits the reverse inequality. The author, attempting to derive the ten necessary stability 
criteria, relied heavily upon the inequality (2.6), which is manifestly in contradiction with 



the second law of thermodynamics. The ten stability criteria, taking t 



of the type d t A < d t B + d t C (see Eqs (3.1), (3.2) and (3.3) in Ref. [l]), have been deduced 
from Eqs (1.5) and from the inequality (2.6). Hence, they are also incorrect. The application 
of the the ten stability criteria to the thermonuclear reactor IGNITOR is also incorrect as 



re form of inequalities 



the discussion revolves around the (incorrect) Eq. (3.1) in Ref. [l]. 
The present comments are organized as follows. The expression for the second order differ- 
ential of entropy, obtained following the original demonstration of Glansdorff-Prigogine is 
illustrated in Section (JTT]) . The link between the Le Chatelier-Braun principle and the local 
stability condition is also shown in this section. The main conclusion of our analysis can be 
found in Section 

II. SECOND ORDER DIFFERENTIAL OF ENTROPY 



Let us introduce the Glansdorff-Prigogine increment 6 Our aim is to derive a set of 
relations coming from the second order quantity 5 2 s. The Gibbs relation for 5s reads: 

T5s = 5u + p5v — PiSNi (3) 

i 

where s indicates the total entropy of the system per unit mass and u and iVj denote 
the energy density per unit mass and the mass fractions, respectively, p, v and are the 
pressure, the specific volume {y = p~ x , with p indicating the total density), and the chemical 
potentials per unit mass, respectively. By simple differentiation of Eq. ([3]) and by expressing 
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dfii in terms of the variables T, p and N i} we obtain the second order quantity 8 2 s [3] 



8 2 s = - 



T 



X 



(Sv) Ni = (--) ST + 

KdTJpM \dp/T,Ni 



ij 
dp 



, vaiv 7 -/T,p,(7v I ) 



where 



(4) 



Here, x is the thermal dilatation coefficient at constant pressure. The index (iVj) means that 
all mass fractions, except N i: are maintained constant. An identical calculation, in which 
the operator 5 is replaced by the time partial derivative dt, yields the equality [?J 



8 t T-%(pu) - p Y J 9t(lnT- l )d t N, + pT~%p d,v - hd t T- l O t p 

\oT / v,Ni Op J T,Ni I \ OT / p,Ni v 

V dNj J T,p,(Ni) 



(5) 



■dv\ T 2 
>Op/T,Ni 



where h stands for the enthalpy per unit mass: h = u + pp~ x . It should be kept in mind that 
in Eq. (jSJ), the independent variables (u,v,Ni) characterize the local state of a dissipative 
system i.e., the convective effects are neglected. According to the thermodynamic stability 
theory, a state is defined to be stable if no evolution starting from the unperturbed state 
can satisfy the requirements of the second law. In the presence of hydrodynamic effects, the 
generalized sufficient condition of local stability, valid for convective as well as for dissipative 

n n 

processes, takes the form [6(, [7[ 



d f S 2 z > ; S 2 z < 



(6) 



where 

S 2 z = 5 2 s - T-^v) 2 = 8T~ x 8u + 5{pT- l )5v - J^^T" 1 )^ - T-\8v) 2 < (7) 

i 

Hence, in case of time-dependent convection processes, a supplementary contribution equal 
to — T _1 (5v) 2 should be added to the second variation of entropy. Of course, in presence of 
hydrodynamic effects the fields, like the temperature, are linked to the velocity through the 
balance equations for mass, energy and momentum. This extra term has been omitted in 
Eq. (1.5) and in Eq. (1.4) of Refs. Jl| and {2 1 respectively. Notice that, as clearly explained 
by P. Glansdorff and I. Prigogine (see Ref. 7J, chapter VI, §7), the terms — T^ 1 (5\) 2 is of 
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the same order as the other contributions and then, by no means can this term be neglected. 
We conclude this comments by recalling that the Le Chatelier-Braun Principle affirms that 
as a result of a variation in one of the factors governing the thermodynamic equilibrium of 
a system, the system tends to adjust to a new equilibrium state counteracting the imposed 
change. It has been proven that, if using 5 2 z as the Lyapunov function the system is locally 
stable then the Le Chatelier-Braun Principle is automatically satisfied j^y]. Using Eq. (jlj) 
with the additional term — T~ 1 (5v) 2 , the second global stability condition reads 

ij 3 i 

where Vj denotes the ith component of the velocity of matter. Let us now consider sys- 
tems subject to time-independent boundary conditions and, in particular, to the additional 
boundary conditions [<9t(vj)]s = 0. Using these boundary conditions, and by taking into 
account the balance equations for mass, energy and momentum, we derive the celebrated 
Universal Criterion of Evolution (UCE) by integrating Eq. (J7J) over the volume occupied by 



the system 



6] 



[y^J k d t X K dV<0 (9) 



where J K and X K indicate the thermodynamic flows and the thermodynamic forces, respec- 
tively. As an example of application, we may re-consider the stability of heat conduction 
in a solid where the temperature is the unique variable. For time-independent boundary 
conditions, it is^easy to check that, by applying the divergence theorem, we obtain again 
inequality (J2]) 



III. CONCLUSIONS 

The mistake in the work published in Refs is the result of a wrong interpretation 

of the equation for the second order differential of entropy. In the case of dissipative fluids 
and plasmas 

A. By no means can the increment 5 be replaced by the substantial derivative d t in the 

equation for 8 2 s; 

B. The inequality for the necessary condition for stability should include the additional 

term — T _1 (5v) 2 . 
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Q and Q 



In the papers 111 and |2| we found the (fundamental) mistake reported in A. and the extra 
term,mentioned in B., is totally absent. In addition Eq. (2.6), deduced by the author in 
Ref . [2] , and largely used in his demonstration, is in contradiction with the general expression 
of the UCE. Since Eqs (1.5) and (2.6), in Ref. 1], are the starting point of the author's 
attempt to derive ten necessary conditions for the stability of non-equilibrium magneto- 
hydrodynamic systems, the final inequalities (3.1), (3.2) and (3.3) (in Ref. jl|), as well as 
the derived criteria for stability, are also incorrect. 

The paper cited in Ref. [3j is an application of the (incorrect) variational results published 

n n 

in Refs |1| and |2J to the thermonuclear reactor IGNITOR. Indeed, the discussion revolves 
around Eq.(l) reported in Ref. jjj [which corresponds to Eq. (3.1) in Ref. jl|]. However, 
as stated above, this equation is incorrect and, as a consequence, the conclusions are also 
incorrect. Anyway, calculations for IGNITOR-plasmas, in the collisional transport regimes, 



have been performed recent 
statements reported in Ref. 



y 



9]. The results are not in contradiction with the general 



10| and with the UCE. 

We conclude with a quotation from the original work of Glansdorff and Prigogine on the 
Universal Criterion of Evolution. 

" The sign, which corresponds to the exact differential of the total change of the entropy 
production, is by no means prescribed by the Universal Criterion of Evolution" . 
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